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ABSTRACT 

Cataclysmic variables (CVs) are binaries in which a compact white dwarf accretes ma- 
terial from a low-mass companion star. The discovery of two planets in orbit around the CV 
HU Aquarii opens unusual opportunities for understanding the formation and evolution of this 
system. In particular the orbital parameters of the planets constrains the past and enables us to 
reconstruct the evolution of the system through the common-envelope phase. During this dra- 
matic event the entire hydrogen envelope of the primary star is ejected, passing the two planets 
on the way. The observed eccentricities and orbital separations of the planets in HU Aqr enable 
us to limit the common-envelope parameter a A = 0.45±0.17or7 = 1.77±0. 02 and measure 
the rate at which the common envelope is ejected, which turns out to be copious. The mass in 
the common envelope is ejected from the binary system at a rate of m = 1.9 ± O.3M0/yr. 
The reconstruction of the initial conditions for HU Aqr indicates that the primary star had a 
mass of A/zams = 1-6 ± 0.2M Q and amzAMS = 0.47 ± 0.04 M Q companion in a a = 25- 
160 R.0 (best value a = 97 R©) binary. The two planets were born with an orbital separation 
of a a = 541 ± 44 R and at, = 750 ± 72 R respectively. After the common envelope, 
the primary star turns into a 0.52 ± 0.01 M helium white dwarf, which subsequently ac- 
creted ~ 0.30 M from its Roche-lobe filling companion star, grinding it down to its current 
observed mass of 0.18 M . 
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1 INTRODUCTION 

The cataclysmic variable HU Aqr currently consists of a 0.80 Mg 
white dwarf that accretes from a 0.18 Mq main-sequence com- 
panion star jSchwopeetai]|201lh . The transfer of mass in the 
tight a — 0.82 R0 orbit is mediate by the emission of grav- 
itati onal waves and the str ong magnetic field o f the accreting 
star 1 Verbunt & Zwaanlll98lh . Since its discovery dSchwope et all 
1 1993 ). irregularities of the observed-calculated variations have 
led to a range of explanations, including the presence of 
circum-binary planets 1 Schwarz et alj 1 20091 : iHorner et all l2012bl : 
iGozdziewski et all 120120 . Detailed timing analysis has eventu- 
ally led_to_the conclusion that the CV is orbited by two plan- 
ets fainse et alj|2012h . a 5.7M Jup planet in a ~ 1205 R© orbit 
with an eccentricity of e = 0.20 and a somewhat more massive 
(7.6 Mj up ) planet in a w ider 1785 R© and eccentric e = 0.38 
orbit ( IHorner et alj|2012bh . Although, the two-planet configuration 
turned out to be dynamicall y unstable on a 1000 — 10,000 year time 
scale (Horner etal.ll2012bl . see also §[4]l, a small fraction of the 
numerical s imulations exhibit lo ng term dynamical stability (for 
model B2 in Wittenm ver et al.l2012l . see Tab.[T]for the parameters). 

It is peculiar to find a planet orbiting a binary, in particu- 



lar around a CV. While pla nets may be a natural consequenc e 
of the formation of binaries dPelupessv & Portegies Zwartll2012l) . 
planetary systems orbiting CVs could also be quite common. 
In part icular because of recen t ly timing re s idual i n NN Ser- 
pentis (Beuermann e t al.1 120101 : IHorner et alj l2012a|). DP Leo- 
nis feeuermann et alj|201lh and OS Virgo dAlmeida & lab lonski 
1201 ll were also interpreted is being caused by circum-CV planets. 

Although the verdict on the planets around HU Aqr (and the 
oth er CVs) remains debated (Tom Marsh private communication, 
and lOian et al. 2010), we here demonstrate how a planet in orbit 
around a CV, and in particular two planets, can constrain the CV 
evolution and be used to reconstruct the history of the inner binary. 
We will use the planets to perform a precision reconstruction of the 
binary history, and for the remaining paper we assume the planets 
to be real. 

Because of their catastrophic evolutionary history, CVs seem 
to be the last place to find planets. The original binary lost proba- 
bly more than half its mass in the common-envelope phase, which 
causes the reduction of the binary separation by more than an or- 
der of magnitude. It is hard to imagine how a planet (let alone two) 
can survive such turbulent past, but it could be a rather natural con- 
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sequence of the evolution of CVs, and its survival offers unique 
diagnostics to constrain the origin and the evolution of the system. 



2 THE EVOLUTION OF A CV WITH PLANETS 

After the birth of the binary, the primary star evolved until it over- 
flowed it Roche lobe, which initiated a common-envelope phase. 
The hydrogen envelope of the primary was ejected quite suddenly 
in this episode dWebbink|[ l984), and the white dwarf still bears the 
imprint of its progenitor: the mass and composition of the white 
dwarf limits the mass and evolutionary phase of its progenitor star 
at the moment of Roche-lobe overflow (RLOF). For an isolated 
binary the degeneracy between the donor mass at the moment of 
RLOF (Mrlof), its radius -Rrlof and the mass of its core M core 
cannot be broken. 

The presence of the inner planet in orbit around HU Aqr 
dHinse et alj|2012l ; iHorner et al.ll2012bl : iGozdziewski et al.ll2012h 
allows us to break this degeneracy and derive the rate of mass loss 
in the common-envelope phase. The outer planet allows us to vali- 
date this calculation and in addition to determine the conditions un- 
der which the CV was bom. The requirement that the initial binary 
must have been dynamically stable further constrains the masses of 
the two stars and their orbital separation. 

2.1 Pre common envelope evolution 

During the CV phase littl e mass is lost from the binary system 
M cv ~ constant (but see ISchenker&Kin3l2002h . and the cur- 
rent total binary mass (A/ cv = 0.9 8 Mr) was not affected by 
the past (and current) CV evolution (Ritterl l2~010l) . The observed 
white dwarf mass then provides an upper limit to the mass of the 
core of the primary star at the moment of Roche-lobe contact, and 
therefore also provides a minimum to the companion mass via 
Mcv - M coro . 

With the mass of the companion not being affected by the 
common envelope phase, we constrain the orbital parameters at the 
moment of RLOF by calculating stellar evolution tracks to mea- 
sure the core mass Af corc and the corresponding radius R(M COIB ) 
for stars with zero-age main-sequence mass Mzams- In Fig.Q] we 
present the evolution of the radius of a 3 Mq star as a function of 
A/c 0r e, which is a measure of time. 

We adopted the Henyey stellar evolution code MESA 
dPaxtonetai]|201lh to calculate evolutiona ry track of stars from 
Mzams = 1 to 8M using AMUSE3 iPortegies Zwart et all 
2009, l2012h to run MESA and determine the mass of the stel- 
lar core. The latter is measured by searching for the mass-shell in 
the stella r evolution code for w hich the relative Hydrogen fraction 
< 1.0" 9 dTauris&Dewil200lh . 

At the moment of RLOF the core mass is A/ corc and the stel- 
lar ra dius J?rlof = -R(Af corc ). Via the relation for the Roche ra- 
dius (Eggleton 1983), we can now calculate the orbital separation 
at the moment of RLOF hrlof as a function of A/rlof - This sep- 
aration is slightly larger than the initial (zero-age) binary separa- 
tion azAMS due to the mass lost by the primary star since its birth 

The Astrophysics Multipurpose Software Environment, or AMUSE, is 
a component library with a homogeneous interface structure, and can be 
downloaded for free at amusecode . org. All the source codes and scripts 
for reproducing the calculations described in this paper are available via the 
AMUSE website. 



A/rlof — Mzams - The long (main-sequence) time scale in which 
this mass is lost guarantees an adiabatic response to the orbital sep- 
aration, i.e. aAItot = constant. 

For each Mzams we now have a range of possible solutions 
for hrlof as a function of M corc and m comp = M cv - M corB . 
This reflects the assumption that the total mass (m comp + M cor c = 
0.98 Mq) in the observed binary with mass M cv is conserved 
throughout the evolution of the CV. In Fig.Q]we present the corre- 
sponding stellar radius .Rrlof and orlof as a function of M core 
for Mzams = 3Mq. This curve for orlof is interrupted when 
RLOF would already have been initiated earlier for that particu- 
lar orbital separation. We calculate this curve by first measuring 
the size of the donor for core mass M coro , and assuming that the 
primary fills its Roche-lobe we calculate the orbital separation at 
which this happens. 



2.2 The common envelope evolution 

During the common envelope phase the primary's mantle is blown 
away beyond the orbit of the planets. The latter responds to this by 
migrating from the orbits in which they were born (semi-major axis 
a a and eccentricity e a , the subscript 'a' indicates the inner planet, 
we adopt a 'b' to indicate the outer planet) to the currently observed 
orbits. Using first order analysis we recognize two regimes of mass 
loss: fast and slow. In the latter case the orbit expands adiabatically 
without affecting the eccentricity: The minimum possible expan- 
sion of the planet's orbit is achieved when the common envelope is 
lost adiabatically. Fast mass loss leads to an increase i n the eccen- 
tricity as well and may even cause the planet to escape dHillsl 19831 : 
IPiiloo et al.ll2012l) . 

A planet born at the shortest possible or bital separation to be 
dynam ically stable will have a a ~ 3ozams dMardling & Aarsethl 
1200 ll) . which is slightly smal ler than the distance at which circum 
binary planets tend to form dPelupessv & P ortegies Z^vartll2012l) . 
In Fig.Q]we present a minimum to the semi-major axis for a planet 
that was born at a a = 3azAMS and migrated by the adiabatic loss 
of the hydrogen envelope from the primary star in the common- 
envelope phase. The planet can have migrated to a wider orbit, but 
not to an orbit smaller than the solid black curve (indicated with 
min(a a )) in Fig.Q] For the 3 Mq star, presented in Fig.Q] RLOF 
can successfully result in the migration of the planet to the observed 
separation in HU Aqr for M core 0.521 Rq, which occurs for 
azAMS HI R©- A core mass Af CO r C > 0.521 Rq would, for a 
3 Mq primary star, result in an orbital separation that exceeds that 
of the inner planet in HU Aqr; in this case the core mass of the 
primary star must have been smaller than 0.521 Mq. 

Another constraint on the initial binary orbit is provided by the 
requirement that the mass transfer in the post common-envelope 
binary should be stable when the companion starts to overfill its 
Roche lobe. To guarantee stable mass transfer we require that 
wicomp ~ M core . The thick part of the red curve in Fig.Q]indicates 
the valid range for the initial orbital separation and core-mass for 
which the observed planet can be explained; the thin parts indicate 
where these criteria fail. 

We repeat the calculation presented in Fig.Q] for a range of 
masses from A/zams = 1 Mq to 8Mq with steps of 0.02 Mq, 
the results are presented as the shaded region in Fig.|2] 
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Figure 1. The orbital separation as a function of core mass for a primary 
star of Mzams = 3Mq. The thin blue curve (bottom) gives the stellar 
radius -Rrlof as a function of its core mass M coro , and the curve di- 
rectly above that (red) gives the initial orbital separation flzAMS f° r which 
Roche-lobe overflow of that binary occurs. Here we adopted a companion 
mass of m C omp = O.98M0 — M core as in HU Aqr. The interruptions 
in the curves indicate the core masses for which Roche-lobe overflow can- 
not occur, because it would already have occurred in an earlier stage of the 
evolution, i.e., at a smaller core mass. The solid black curve (top) gives 
the minimal orbital separation of a planet born at a a = 3azAMS f° r with 
the orbit was adiabatically expanded due to the mass loss in the common 
envelope A^zAMS — M COTC . The horizontal dotted curve gives the sepa- 
ration at which the inner planet around HU Aqr was observed. For viable 
solutions the solid black curve should remain below the horizontal dotted 
curve. The thick parts of the red curve indicate where the zero-age binary 
complies to the most favorable conditions for engaging RLOF, surviving 
the common-envelope and produce a planet that can migrate to at least the 
observed separation for the inner planet in HU Aqr. 



3 EFFECT OF THE COMMON ENVELOPE PHASE ON 
THE PLANETARY SYSTEM 

The response of the orbit of the planet to the mass loss depends on 
the total amount of mass lost in the common envelope and the rate 
at which it is lost. Numerical common-envelope studies indicate 
that f or an in-spiraling binary m ~ 2.OM0/yr l lRicker & Taaml 
120121) . At this rate the entire envelope Mrlof - M COTe ~ 0.46- 
5.8 Mq is expelled well within one orbital period of the inner 
planet, which leads to an impulsive response and the possible loss 



(for A/rlof - M c , 



1.92 Mq) of the planet. The fact that the 



HU Aqr is orbited by a planet indicates that at the distance of the 
planet m a < (A/rlof - M core )/P p ianct ~ lM Q /yr. The ec- 
centricity of the inner planet in HU Aqr (see Tab.[TJ can be used 
to further constrain the rate at which the common-envelope was 
lost from the planetary orbit. The higher eccentricity of the outer 
planet indicates a more impulsive response, which is a natural con- 
sequence of its wider orbits with the same m. This regime between 
adiabatic and impulsive mass loss is hard to study analytically (Li 
l2008h . 



3.1 The response of the inner planet 

We calculate the effect of the mass loss on the orbital parameters 
by numerically integrating the planet orbit. The calculations are 
started by selecting initial conditions for the zero-age binary HU 
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Figure 2. Distribution of initial conditions azAMS an d A-izAMS. and the 
resulting core mass M COTC (bottom panel) that successfully reproduce the 
CV HU Aqr. In the top panel we also present the results of our analysis 
for the two planets. The circles give the initial semi-major axis of the inner 
planet. The triangles give the initial semi-major axis for the outer planet. 
The symbols are colored (red for the inner planet and blue for the outer) if at 
least 10 out of 20 calculations for the orbital integration over 1 Myr turn out 
to be dynamically stable (see|4). With these initial conditions both planets 
migrate to within 1 % of their observed orbital period with an eccentricity of 
e a = 0.20 ±0.01 ande 6 = 0.33 ± 0.03 with m a = 0.15 ± 0.01 M /yr. 



Aqr -Mzams , Izams an d consequently Af corc - from the available 
parameter space (shaded area) in Fig.[2] and integrate the equations 
of motion of the inner planet with time. Planets ware assumed to 
be born in a circular orbit (e a — 0) in the binary plane with semi- 
major axis a a . 

The equations of motions ar e integrated using th e high- 
order symplectic integrator Huayno dPelupessv et al.ll2012f) via the 
AMUSE framework. During the integration we adopt a constant 
mass-loss rate m applied at every l/100th of an orbit, and we con- 
tinued the calculation until the entire envelope is lost (see § !2.2l and 
Fig.|2j, at which time we measure the final semi-major axis and ec- 
centricity of the planetary orbit. During the integration we allow 
the energy error to increase up to at most AE/E — 10~ 13 . 

By repeating this calculation while varying a a and rh a we it- 
erate (by bisection) until the result is within 1% of the observed 
aHUAqr(o) and eHUAqr(a) of the inner planet observed in HU Aqr. 
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Table 1. Reconstructed and observationally constrained parameters for HU Aquarii. The parameters at zero-age (column 2 and 3) and directly after the 
common-envelope phase (columns 4 and 5) are derived by means of reconstructing the CV evolution. The best comparison is achieved for a mass loss rate in 
the common-envelope from the binary system of m = 1.9 ± 0.3 Mg/jr (or ih a = 0.15 ± 0.01 Mp/yr f rom the orbi t of the first planet). Parameter s that 
we were unable to constrain are printed in slanted font, observed parameters are from (Schwope et al. 201 1 ; Hinse et al. 2012; Gozdziewski et al. 2012j) and 
are printed in boldface. 



parameter Zero age After CE Today 



M/M 1.6 ±0.2 0.52 ±0.01 0.80 ± 0.04 

m/Mg 0.47 ±0.04 0.47 ± 0.04 0.18 ±0.06 

a/R 25-160 0.867-2.0 0.8 

e 0.0 0.0 0.0 

planet a planet b planet a planet b planet a planet b 

a/R 541 ±44 750 ± 72 1204 ± 12 1785 ± 18 1204 1785 

m/Mj U p 5.7 7 5.7 7 5.7 7 

e 0.0 0.0 0.20 ±0.00 0.33 ± 0.03 0.20 0.38 



The converged results of these simulations are presented in Fig.[2] 
(circles), and these represent the range of consistent values for the 
inner planet's orbital separation a a = 183-752 R as a function 
of AizAMS = 1-8 M and consistently reproduce the observed in- 
ner planet when adopting m a — 0.124-0.267 M /yr. The highest 
value for rh is reached for A/zams = 2.85 M at an initial orbital 
separation of a a = 427 R . 

The orbital solution for the inner planet is insensitive to 
the semi-major axis of the zero-age binary ozams (for a fixed 
A/zams), and each of these solutions were tested for dynamical 
stability, which turned out to be the case irrespective of the initial 
binary semi-major axis (as discussed in §[4j. 

3.2 The response of the outer planet 

We now adopt the in § l3,ll measured value of m to integrate the orbit 
of the outer planet. The effect of the mass outflow on the planet is 
proportional t o the square of the d ensity in the wind at the location 
of the planet dKudritzki et "al] |l992l) . We correct for this effect by 
reducing the mass loss rate in the common envelope that affects the 
outer planet by a factor (a,b/a a ) 3 ^ 2 ■ 

We use the same integrator and assumptions about the initial 
orbits as in § 13.11 but we adopt the value of m from our reconstruc- 
tion of the inner planet (see §[34). To reconstruct the initial orbital 
separation of the outer planet at, we vary this value (by bisection) 
until the final semi-major axis is within 1% of the observed orbit 
(see Tab.[T|l. The results are presented in Fig.[2](triangles). The post 
common-envelope eccentricity of the outer planet then turn out to 
be e b = 0.38 ± 0.07. 



4 STABILITY OF THE INITIAL SYSTEM 

After having reconstructed the initial conditions of the binary sys- 
tem with its two planets we test its dynamical stability by inte- 
grating th e entire system numeri cally for 1 Myr using the Huayno 
integrator iPelu pessv et al.l2012l) . To test the stability we check the 
semi-major axis and eccentricity of both planets every 100 years. 
If any of these parameters change by a factor of two compared to 
the initial values or if the orbits cross we declare the system un- 
stable, otherwise they are considered stable. The calculations are 
repe ated with the 4th orde r Hermite predictor-corrector integrator 
ph4 dMcMillan et al.l2012l) within AMUSE to verify that the results 
are robust, which turned out to be the case. We then repeated this 
calculation ten times with random inital orbital phases and again 



with a 1% Gaussian variation in the initial planetary semi-major 
axes. In Fig.|2]we present the resulting stable systems by coloring 
them red (circled) and blue (triangles), the unstable systems are 
represented by open symbols. 

From the wide range of possible systems that can produce HU 
Aqr only a small range around A/zams = 1.6 ± 0.2 M turns out 
to be dynamically stable. The eccentricity of the outer orbit of the 
stable systems (which ware stable for initial conditions within 1%) 
et = 0.32 ± 0.02, which is somewhat smaller than the observed 
value for HU Aqr faorner et al]|2012bl . e = 0.38 ± 0.16). These 
values are obtained with rn a — 0.15 ± O.OlM /j/r. The small 
uncertainty in the derived value of m is a direct consequence of its 
sensitivity to e a and the small error on Mzams from the require- 
ment that the initial system is dynamically stable. 



5 DISCUSSION AND CONCLUSIONS 

We have adopted the suggestive results from the timing analysis of 
HU Aqr, that the CV is orbited by two planets, to reconstruct the 
evolution of this complex system. A word of caution is well placed 
in that these observations are not confirmed, and currently under 
debate (Tom Marsh private communication, and comments by the 
referee). However, The predictive power that such an observation 
would entail is interesting. The possibility to reconstruct the initial 
conditions of a CV by measuring the orbital parameters of two cir- 
cum binary planets is a general result that can be applied to other 
binaries. For CVs in particular it enables us to constrain the value of 
fundamental parameters in the common-envelope evolution. This in 
itself makes it interesting to perform this theoretical exercise, irre- 
spective of the uncertainty in the observations. On the other hand, 
the consistency between the observations and the theoretical anal- 
ysis give some trust to the correctness of these observations. 

The presence of one planet in an eccentric orbit around a CV 
allow us to calculate the rate at which the common-envelope was 
lost from the inner binary. A single planet provides insufficient in- 
formation to derive the initial mass of the primary star, but allows 
us to derive the initial binary separation and planetary orbital sep- 
aration to within about factor of 5, and the initial rate of mass loss 
from the common envelope to about a factor 2. A second planet 
can be used to further constrain these parameters to a few per cent 
accuracy and allows us to make a precision reconstruction of the 
evolution of the CV. 

We have used the observed two planets in orbit around the CV 
HU Aqr to reconstruct its evolution, to derived its initial conditions 
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(primary mass, secondary mass, orbital separation, and the orbital 
separations of both planets) and to measure the rate of mass lost 
in the common-envelope parameters m. By comparing the binary 
parameters at birth with those after the common-envelope phase we 
subsequently calculate the two parameters aX and 7. 

The measured rate of mass loss for HU Aqr of rh a — 
0.15 ± 0.01 Mq /yr from the inner planetary orbit, which from 
the binary system itself would entail a mass-loss rate of m = 
1.9 ± 0.3 Mq /yr, when we adopt the initial binary to have a semi- 
major axis of azAMS — 97 Rq, which is bracketed by our derived 
range of ozams = 25-160 R©. This is consistent with a mass-loss 
rate of rh ~ 2 M© / y r from numerical common-envelope studies 
dRicker & Taamll2012l) . 

By adopting that the binary survives its common envelope at a 
separation between ~ 0.87 R0 (at which separation the secondary 
star will just fill it's Roche-lobe to the white dwarf) and ~ 2 Rq 
(for gravitational wave radiation to drive the binary into Roche-lobe 
overflow within lOGyr), we derive the value of aX — 0.2-2.0 (for 
azAMS — 97 R© we arrive at aX ~ 0.45 ± 0.17). This value is a 
bit small compared to numerous earlier studies, which tend to sug- 
gest aX ~ 4.0. The alternative 7-formalism for common-envelope 
ejection JNelemans et alj200oT) gives a value of 7 = 1.63-1.80 (for 
azAMS — 97 R© we arrive at 7 ~ 1.77 ± 0.02 ), which is consis- 
tent w ith the determination of 7 in 30 other CVs jNelemans & TouJ 
120051) . 

The inner planet in HU Aqr formed at (3.2-20.6)cizams, with 
a best value of 5.3±0.45cizams, which is consi stent with the plan- 
ets found to orbit other bin aries, like Kepler 1 6 dDovle et alj|201lh 
and for Kepler 34 and 35 dWelsh et alj|2012l) . although these sys- 
tems have lower primary mass and secondary mass stars. 

It seems unlikely that more planets were formed inside the or- 
bit of the inner most planet, even though currently there is sufficient 
parameter space for many more stable planets; in the zero-age bi- 
nary there has not been much room for forming additional planets 
further in. It is however possible that additional planets formed fur- 
ther out and those, we predict, will have even higher eccentricity 
than those already found. 
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